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We consider a Bose-Einstein condensate driven by periodic (5-kicks. In contrast to first-order descriptions, 
which predict rapid, unbounded growth of the noncondensate in resonant parameter regimes, the consistent treat- 
ment of condensate depletion in our fully-time-dependent, second-order description acts to damp this growth, 
leading to oscillations in the (non)condensate population and the coherence of the system. 
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Central to the concept of a weakly interacting atomic Bose- 
Einstein condensate (BEC) is that each of the N component 
atoms can be considered to be in approximately the same mo- 
tional state; this is manifest through the mean-field descrip- 
tion of zero-temperature (T - 0) BEC dynamics with the 
Gross-Pitaevskii equation (GPE), which takes the form of a 
single-particle Schrodinger wave equation with an additional 
cubic nonlinear term (TJ. At finite temperature there is ther- 
mal depletion of the condensate, which can be theoretically 
accounted for in a variety of ways [2]. Even at T — in a 
system of finite size there is always a finite noncondensate 
fraction |[T] |3), and one expects dynamics within the BEC 
to cause significant particle transfer from the condensate to 
the noncondensate fraction under quite general circumstances 
When rapid, such dynamical depletion has commonly 
been supposed to presage destruction of the BEC as a coher- 
ent entity, however previous studies have been hampered by 
the absence of a self-consistent treatment 13. Applying 
the number-conserving approach of Gardiner and Morgan J8), 
fully dynamically to second order, we have carried out the first 
such self-consistent treatment of a specific example system. 

Our chosen test-system is based on the quantum 5-kicked 
rotor ll9 UT2l . a paradigm quantum-chaotic system in which 
periodic driving leads to complex behavior, including dynam- 
ical localization J9] [12l and quantum resonances (associated 
with ballistic increase in the kinetic energy ll9 fflTII ). Atom- 
optical realizations of such systems lfTTI - [T3l comprise an ex- 
citing area of research into quantum-chaotic phenomena; ex- 
tension into the regime of BECs has also become an active 
area of research, in which several new phenomena have been 
predicted JS] [7] [14UT61 . In the mean-field approximation, 
the GPE nonlinearity can strongly influence 5-kicked-rotor- 
BEC dynamics [14|; in particular, the structure of quantum 
resonances was recently elucidated [16|, revealing previously 
unobserved resonance profiles with a sharp asymmetric cut- 
off. The noncondensate fraction can be accounted for using 
a number-conserving approach (3] |4] [8] [T7) : such approaches 
contain a fixed atom number (appropriate to atomic BEC ex- 
periments ID), and nonlocal terms [/ and Q in Eqs. ^ and 
Q] which, within the second-order, linear-response treatment 
used by Morgan et al. Ifl8l . have proved vital in explaining 
the temperature dependence of collective excitation spectra 




FIG. 1. (Color online) (a) (5-kicked-rotor-BEC: (i) toroidally-trapped 
BEC; (ii) <5-kicking potential, (b) Evolution of condensate and non- 
condensate fractions n c = N C JN and n, = N,/N in (i) first- and 
(ii) second-order number-conserving descriptions (N = 10 4 , gr = 
2.5 x 10~ 4 , T p = 9.255, k = 0.5). In (iii) we show the coherence 
measure C [Eq. j7]l] in the second-order description, and the fidelity, 
F [Eq. {8}], of the condensate mode between descriptions. 



observed at JILA and MIT Ifl9l . The first-order, number- 
conserving description of Castin and Dum [4] [consisting of 
modified Bogoliubov-de Gennes equations (MBdGE) coupled 
to the GPE] has been the approach of choice in (5-kicked BEC 
systems 017]. This has revealed a general tendency towards 
rapid growth in the number of noncondensate atoms N t . Such 
behavior is directly linked to linear dynamical instabilities 
(i.e., sensitivity to initial conditions in the linearized regime) 
in the GPE, enabled by its nonlinearity [5]; the presence of 
such instabilities is a generic feature of most nonlinear sys- 
tems. Conversely, growth in N, should match depletion in the 
condensate number N c , and as atoms transfer from N c to N,, 
qualitatively one expects mean-field interactions and hence 
further transfer to "switch off" at some stage. Whether this 
occurs before the destruction of the condensate has remained 
an open question, as a linearized, first-order description treats 
the condensate as an effectively undepletable "particle bath," 
feeling no effect from the noncondensate. 

In this Letter we answer this question in the affirmative, 
using the first fully time-dependent application of the second- 
order, number-conserving description of Gardiner and Mor- 
gan OD, for a 5-kicked-rotor-BEC in a quasi-one-dimensional 
ring trap [Fig.[T|a)], an ideal test system to explore generic is- 
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sues of dynamically induced condensate depletion. This self- 
consistent description consists of a generalized GPE (GGPE) 
coupled to the MBdGE; it explicitly preserves N — N c + N, 
and the orthogonality of the condensate and noncondensate, 
includes mutual interactions, and allows free transfer of popu- 
lation between the two. We solve these equations numerically 
and explore resonant parameter regimes which, to first order, 
lead to rapid, unbounded growth of the noncondensate. Our 
principal finding is the damping of this growth in the second- 
order description [Fig.[TJb)(i,ii)]. We also compute the coher- 
ence of the system and the departure of the second-order de- 
scription from the GPE [Fig.[TJb)(iii)] for varying N. We show 
that, despite considerable differences in dynamics between the 
descriptions around resonant parameter regimes, the GPE ac- 
curately predicts the location of these parameter regimes; the 
cut-off identified in [16| is qualitatively preserved, however 
the accompanying exponential oscillations are strongly modi- 
fied for experimentally realistic atom numbers. 

We consider N bosonic atoms of mass M, held in a toroidal 
potential V T (p,z) = Mcu 2 [(p - R) 2 + z 2 ]/2 [Fig. |jja)(i)], 
interacting with s-wave contact interactions, and subject to 
a temporally and spatially periodic driving potential V(0, t) 
[Fig.[TJa)(ii)]; toroidal potentials similar to Vj can be created 
and precisely controlled using all-optical methods ll20l . As- 
suming sufficiently strong radial and axial confinement, and 
harmonic length r = y/h/Mu> <*: R CD . we reduce the sys- 
tem Hamiltonian to a dimensionless (length unit R, time unit 
MR 1 IK), one-dimensional form [16|: 



+ V(8, t) + —^(8)^(8) 
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where the field operators obey bosonic commutation rela- 
tions [¥(0), 4^(0')] = 6(8 - 8'). The interaction strength 
gr = 2a s R/r 2 , where a s is the s-wave scattering length. As 
in ll7l [T6ll . we model the driving potential as a train of 5-kicks, 
V(0, t) - kcos(0) 2^=0 _ n Tp), with (dimensionless) kick- 
ing period T p . This may be approximated in experiment using 
short pulses of off -resonant laser light lfT0lfT2ll22ll . 

We define the condensate mode if/(8) (with creation oper- 
ator a c ) as the eigenfunction of the single-body density ma- 
trix ($r(ffy¥(Q)) with the largest eigenvalue N c (the num- 
ber of condensate atoms), to which it is normalized, i.e., 
f d8\i//(8)\ 2 = N c = (ala c ). Hence, we expand the field 
operator as ¥(0) = a c tff(6)/ a/A^ + &¥(ff), where &¥(6) = 
J d8Q(6,ffyV{ff) describes the field component orthogonal 
to the condensate, and the projector Q(8, 0') = 6(8 - 8') - 
^(8)^*(8')IN C . We propagate if/(8) with the GGPE: 
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where A 2 is a real-valued scalar constant with a role similar to 
the chemical potential in grand-canonical treatments 11231 . and 
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Introducing the number-conserving fluctuation operator 
A(0) = al8%(8)l EJ, we define n(8,8') = <A t (0')A(0)) 
and m(8, 0') = (A(0')A(0)) (the noncondensate normal and 
anomalous averages), and 

f(0)=^- J d8'W8')\ 2 [H8,8')m) + r(0')m(O',O)], (4) 

which ensures orthogonality of the condensate from the non- 
condensate component. The dynamics of tfr(8) are therefore 
coupled to those of A(0), A t (0). We decompose A(0) = 
^dbkum+b-kU-m+^vli^+bl^m, where b k , b\ are 
bosonic quasiparticle operators, and the quasiparticle modes 
are normalized to J dd[\uk(8)\ 2 - |v'i;(0)| 2 ] = 1, and choose 
all time-dependence to be within the quasiparticle mode func- 
tions. These are then propagated by the MBdGE: 
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where L(0,0') = 6(8 - 8')[H G? (8') - Aq] + 

g T fd8"Q(8,8")M0")\ 2 Q(8",0'), M(8,8') 

g T Jd8"Q(8,8")ip(8") 2 Q*(8",8'), and A [the system 

ground state value of (1/AQ J ' dd\p*(8)H GV (8)if/(8)] is a 
lower-order approximant to Az- At T = 0, we can thus at 
all times express n(0,0') = 2X1 [v k {8)v\{8') + V- k (0)v*_ k (8')l 
m(8,8') = E£ 1 [a t (fl)vj(0') + M-^vl^)]. The k index 
quantifies the momentum associated with the equilibrium 
quasiparticle eigenmodes, prior to application of V(0, t). 
The coupling of Eq. (|2]) and Eq. |5]l constitutes the second- 
order, number-conserving description of Gardiner and 
Morgan |8], and represents the equation of motion for 
( P(0) expanded to second order in the fluctuation operator 
A(0) OS)- To first-order @ [T7) one obtains the GPE 
idifr(8)/dt = [H GP (8) - Ao]ifr(8), coupled to the MBdGE 
[Eq. (BJ], as used in previous time-dependent studies of 
noncondensate dynamics of (5-kicked BECs 00. The GPE 
alone constitutes a zeroth-order description, although it may 
be possible to infer higher-order processes from a pure GPE 
treatment |[T6l : unlike GPE plus MBdGE, this is at least an 
internally consistent theoretical description |2]|8). 

We take the T = equilibrium state (without driving) as our 
initial condition. The initial condensate mode is therefore spa- 
tially homogeneous: iff — y]N c l2n. This sets Aq — gjN c !2n, 
and the initial stationary quasiparticle modes 
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where A k = A_ A . = (1 + 4A /k 2 )- 114 . Hence, N t = N - 
N c = j d8h{8,8) = (l/2)2™ =1 (Ai - A" 1 ) 2 , and we set A 2 = 
(gT/27i)[N - 1 + 2itli(A 2 - 1)]. To numerically determine a 
self-consistent T — solution to Eqs. (|2]i and PJ, for given 
values of and gj, we set N c = N, and then; (a) calculate Ak 
up to a cut-off momentum \k\ = m; (b) determine N, from the 
Ak', (c) make the replacement N c = N - N,. We repeat steps 
(a)-(c) until convergence. To determine the driven dynamics, 
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FIG. 2. (Color online) BEC response in the second-order description 
(A' = 10 4 , k = 0.5): (a) relative population of k t momentum 
modes among all atoms, (1 - no), and noncondensate fraction, (n t ); 
(b) coherence measure (C) [Eq. |7)]; and (c) (1 - n ) as predicted by 
the GPE, minus its value in the second-order description. Averages 
are taken over the first 100 kicks. 
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FIG. 3. (Color online) Comparison of first- and second-order de- 
scriptions close to a nonlinear resonance explored in [16| [grN = 
2.5, T p = 6.12, k = 0.5: circle in Fig.[2ja)]. Condensate and noncon- 
densate fractions n c and n t , coherence measure C [Eq. R)], and the 
fidelity of the condensate mode between descriptions, F [Eq. fS)], 
are shown. Columns correspond (left to right) to N = 10 4 , 10 s , and 
10 12 ; agreement of the initial growth in n, between the second-order 
(dependent on g T and N) and first-order (dependent on grN) descrip- 
tions over such a range is a useful test of the second-order numerics. 



we use a Fourier pseudospectral split-step method 11251 : all 
simulations are converged in timestep, grid size, and quasi- 
particle cut-off momentum m. 

In Fig.[TJb) we plot the condensate and noncondensate frac- 
tions (n c = N c /N and n, = N,/N) for parameters which, in the 
first-order description, lead to rapid growth of the nonconden- 
sate (becoming unphysical after ~ 20 kicks) [Fig. |TJb)(i)] . In 
the second-order description [Fig. |TJb)(ii)] the "back-action" 
of the noncondensate rapidly damps out this growth, leading 
instead to complementary oscillations in n, and n c . We also 
track the overall coherence of the system through 

c = Jj ded&gi{e,er)gi{&,G), (7) 

where g x {6, 8') = (6)) IN is the first-order correlation 

function, and compare the evolution of if/ in the GGPE with 
the GPE prediction (iAgpe) through the fidelity 

F= \fder GPE m(Q)\ 2 

NN C 

The quantity C equals unity only in the limit of a pure conden- 
sate, where the noncondensate fraction is exactly zero (i.e., 
the single-body density matrix is exactly factorizable). The 
GGPE then reduces to the GPE, and F = 1 . Otherwise both 
C and F take values between zero and unity. In Fig. [TJb)(iii) 
we observe that C follows the oscillations of n, closely, while 
F shows larger amplitude oscillations with revivals. Simi- 
lar behavior persists across the T p —gr parameter space: in 
Fig. |2|a) we show the time averaged response to weak driv- 
ing (k = 0.5), by plotting 1 - «o averaged over the first 100 
kicks; here = (aiak)IN, where SI creates an atom with 
momentum k. The structure of this response over the range 
of Fig. Sa), modeled with the GPE alone, was recently elu- 
cidated [16|: the response is dominated by linear resonances 



corresponding to the first two primary quantum resonances of 
the 5-kicked rotor as Sr — > BlOi Higher-order linear res- 
onances generally decay with increasing g T , while nonlinear 
resonances appear, with no analog in the linear regime fl6l . 
To first-order (GPE plus MBdGE) all these resonant areas of 
parameter space are associated with rapid growth of the non- 
condensate fraction n, due to linear instabilities in the GPE 
dynamics Q. In contrast, we find that in the second-order 
description (GGPE plus MBdGE) this growth is damped out: 
throughout Fig. |2ja) the 100-kick average of n, remains be- 
low 0.6. However, (n t ) is still strongly enhanced in parameter 
regimes with significant resonant response [large (1 - rcoX as 
shown in Fig.[2|a)]. The predictions of the second-order de- 
scription then differ considerably from the standalone GPE 
description, as shown by the behavior of (C) [Fig. |2|b)], and 
the difference in response (as measured by (1 - no)) between 
descriptions [Fig. [2jc)]. Nonetheless, we find that the reso- 
nances are located in the same regions of parameter space in 
both descriptions, and that the asymmetric profiles and sharp 
cut-offs seen in [ 16] remain. Away from these resonances the 
GPE agrees well with the second-order description, up to the 
100 kicks we consider. 

In Figs.[3]and|4]we compare the first- and second-order de- 
scriptions, for varying N but fixed gjN, close to a nonlinear 
resonance studied in [16| [circle in Fig.^a)]. Figure [3] shows 
that the dynamics in the second-order description match the 
GPE for times which increase with N. This increase is slow, 
however; for realistic N (<s 10 8 ) the loss of coherence, unac- 
counted for in the GPE and measured by decay in C, quickly 
becomes significant. Furthermore, in Fig. |4]we see that, on 
the same timescales associated with significant decay in C, 
the dynamics of the relative populations n^, as studied in ifTBI 
using the GPE, noticeably differ in our second-order descrip- 
tion. Compared to the first-order description, rapid growth 
of the noncondensate begins at the same time in our second- 
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FIG. 4. (Color online) Relative populations, n k , of low momen- 
tum modes in the first- and second-order descriptions (parameters 
as in Fig. |3j. In the second-order description n k is shown among 
all atoms and among the noncondensate atoms; in the first-order de- 
scription n k is shown among condensate atoms (as a bar chart) and 
among the noncondensate atoms (with a three-period moving aver- 
age). Columns correspond (left to right) to N = 10 4 , 10 s , and 10 12 . 

order description. However, transfer of population to the non- 
condensate is driven by, and sensitive to, atom-atom interac- 
tions. Hence, decreasing population of the condensate, consis- 
tently accounted for in the second-order description, reduces 
the mean-field interactions, and hence the rate of population 
transfer. We observe population oscillations between con- 
densate and noncondensate fractions, accompanied by oscil- 
lations in the coherence C and fidelity F [Pig. [3) . In Fig.|4]we 
also observe the exponential oscillations in «2 + «-2 reported 
in 1 16 1; however, for realistic atom numbers the frequency of 
these oscillations is quickly increased by the presence of a 
significant noncondensate fraction. 

We have applied a fully dynamical, second-order, number- 
conserving approach |8| to the 5-kicked-rotor-BEC. In con- 
trast to a first-order approach, we observe that rapid growth of 
the noncondensate in resonant parameter regimes is damped 
by our consistent treatment of the condensate population and 
condensate-noncondensate interactions. Although our de- 
scription leads to different dynamics around resonant parame- 
ter regimes, these regimes occur where the GPE predicts them. 
Furthermore, our description retains the cut-offs, but will typi- 
cally strongly modify the exponential oscillations predicted by 
the GPE. Extension of our second-order description to other 
BEC systems offers exciting possibilities for future research. 
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